Charge carrier induced barrier height reduction at organic heteroj unctions 
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In order to provide an accurate theoretical description of current density voltage (J — V) characteristics of 
an organic heterojunction device over a wide range of electric fields at various temperatures, it is proposed that 
an accumulation of charge carriers at the heterojunction will lead to a reduction in the barrier height across the 
heterojunction. Two well-known hole transporting materials, 4,4,4-Tris(N-3-methylphenyl-N-phenyl-amino) 
triphenylamine (MTDATA) and N,N-diphenyl-N,N-bis(l-naphthyl)(l,l-biphenyl)-4,4diamine (NPB) were used 
to fabricate unipolar heterojunction devices. It is found that the J — V characteristics depends strongly on applied 
bias. The simulated J — V characteristics of the heterojunction device, with the modified injection model, are 
found to be in excellent agreement with the experimental data. 

PACS numbers: 71.10.-w, 71.23.-k, 72.20.Ee, 77.80.Le, 73.20.-r, 77.84Jd 



Heterojunctions formed by two semiconductors has been 
the foundational technology for many modern electronic de- 
vices since 1950s [1]. The first high efficiency organic light- 
emitting-diode (OLED) structure based on a stack of multi- 
layer organic materials [2] has become a standard device plat- 
form in commercial production. Recently, 5% power con- 
version efficiency from organic photovoltaic (OPV) cells has 
been demonstrated in small-molecule heterojunction devices 
and in polymeric blend structures [3, 4], and thus low cost 
OPV cells with current printing technology become commer- 
cially attractive. Despite those technological achievements 
with the applications of organic heterojunctions, the studies 
and understanding of the physics of the charge injection pro- 
cess across organic heterojunction are still very limited. 

There are two key factors that control the injection current 
at a heterojunction. One is the energy level alignment which 
determines the barrier height. According to recently reported 
photoemission studies, the interface alignment between two 
undoped organic materials often agrees with the vacuum level 
alignment rule. However, when one side of the heterojunction 
has extrinsic carriers, introduced either by chemical doping or 
electrical doping due to charge carrier accumulation, a consid- 
erable re-alignment of energy levels at the interface has been 
observed [5, 6, 7, 8]. We have also reported that the injection 
current across an organic heterojunction device can be tuned 
by inserting a thin 2 nm chemically-doped interlayer with dif- 
ferent doping concentrations at the interface [9]. The other 
factor is the nature of the charge injection process. It is gener- 
ally believed that the charge transport in organic materials is 
governed by the Miller- Abrahams type hopping process [10] 
between localized transporting sites. Recently, Arkhipov et 
al [11] have proposed a two step hopping model to describe 
the charge carrier transport across an interface of disordered 
organic dielectric. The injection model considers the jump 
rate of a carrier from the center of the Gaussian density-of- 
state (DOS) in the injecting layer into the DOS of the accept- 
ing layer. However, the model does not include the energy 
distribution and occupation of carrier sites at the interface. 
Woudenbergh et al [12] found that measurements of a poly- 



meric heterojunction device do not match the injection model. 
They modified the model to account for the charge filling ef- 
fect of holes in the injecting layer by assuming the carriers are 
injected from the quasi-Fermi level rather than from the center 
of the highest-occupied-molecular orbital (HOMO). 

In this letter, systematic experimental current density volt- 
age (7 — V) characteristics of organic heterojunction devices 
made from relatively well-understood organic molecules were 
studied. In order to describe the experimental results, we 
found that it is essential to include two critical parameters in 
the existing theory: a) energy distribution of carriers at the 
injecting layer interface, and b) dynamic barrier height, (j) v , 
associated with quasi-Fermi level shifting. 

The approach in this experiment is to use a unipolar de- 
vice where the current is mainly determined by the charge 
injection process at the organic heterojunction. This re- 
quires that the contact resistance at the metal-organic (MO) 
interface for charge injection and the number of trap states 
in the bulk organic material are negligible. This can be 
achieved by using 4,4,4-Tris(N-3-methylphenyl-N-phenyl- 
amino) triphenylamine (MTDATA) and N,N-diphenyl-N,N- 
bis(l-naphthyl)(l,l-biphenyl)-4,4diamine (NPB) to fabricate 
the heterojunction device. MTDATA and NPB are proto- 
typical hole-injecting and hole-transporting materials used in 
OLED technology. Due to the difference in HOMO levels 
of MTDATA (5.0-5.1 eV) [13] and NPB (5.4-5.6 eV) [6, 7], 
there exists an energy barrier of around 0.5 eV for hole injec- 
tion at the MTDATA / NPB interface. It has also been demon- 
strated that an Ohmic contact is formed at the indium tin oxide 
(ITO) and MTDATA interface [9, 14]. Moreover, the trap-free 
space-charge-limited conduction behavior in single layer de- 
vices supports that there are a negligible number of trap states 
in both materials [9, 14, 15]. 

A hole-only heterojunction device with a structure of anode 
/ injecting layer / accepting layer / cathode was fabricated for 
this study. Pre-patterned ITO coated glass substrate was used 
as the anode. The substrate was cleaned by a sequence of ul- 
trasonic solvent baths and followed by UV-ozone treatment. 
A 320 nm layer of MTDATA was thermally evaporated on top 
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FIG. 1: (a) The normalized energy distribution of the carrier density, 
p,„,i(£; nt i), at the injecting layer side of the interface, (b) The nor- 
malized fractional injection current densities, J(E[ nt \ ) contributed by 
the carriers at different energy levels. The results are compared with 
different values of the quasi-Fermi level (from —0.2 eV to —0.4 eV). 



of ITO as the injecting layer. It was then covered by 390 nm of 
thermally evaporated NPB as the accepting layer. Finally, sil- 
ver was evaporated on top of NPB to form the cathode. The or- 
ganic and metal layers were fabricated in a single vacuum sys- 
tem with a base pressure of 10~ 8 Torr. On the cathode side, a 
large energy difference between the lowest unoccupied molec- 
ular orbital (LUMO) of NPB (2.5 eV) and the work function 
of Ag (4.3 eV) effectively restricts electron injection into the 
device. The field and temperature dependent hole mobilities 
of individual organic material were characterized by the time- 
of-flight (TOF) technique. Details of TOF measurement has 
been report elsewhere [9]. The measured mobilities were fit- 
ted to the Gaussian-disordered model (GDM) [16] to obtain 
the energetic disorder of HOMO of Omtdata = 93 meV and 
Onpb = 90 meV. 

In order to compare the injection model to the measured 
J — V characteristics of an heterojunction device, the elec- 
tric field at the interface, fj nt , and the charge carrier densi- 
ties at both sides of the interface region, the injecting layer 
side, pimi, and the accepting layer side, pi n a, have to be de- 
termined. The complete two layer organic system including 
contacts was modeled using ATLAS [17], to solve the stan- 
dard electric field, drift, and diffusion equations used in semi- 
conductor modeling. The iteration method was similar to the 
method used in Ref. [12]. The applied voltages and mea- 
sured current densities were used as the input parameters to 
a steady-state model of the device. The field and tempera- 
ture dependent mobilities measured by TOF were used in the 
simulation [9]. The hole-injecting electrode was modeled us- 
ing a Schottky contact with the barrier height chosen to give 
a density of holes equal to 10 21 cm~ 3 at the anode / injecting 
layer interface. The cathode / accepting layer interface was 
specified as Ohmic (i.e. the hole concentration is pinned to 
the equilibrium value at this contact). ATLAS does not con- 
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FIG. 2: A schematic diagrams of energy level alignment at the het- 
erojunction with the boundary condition of (a) vacuum level align- 
ment and (b) Fermi level alignment. The lower graphs are the corre- 
sponding alignment of the Gaussian density-of-state (DOS) at both 
sides of the interface. The shaded region represents those states oc- 
cupied by electrons. 



tain the hopping injection model. Instead, we used a thermal 
injection model at the organic heterojunction and adjusted the 
parameters to give our measured current for each applied volt- 
age. This is equivalent to modeling the two layers separately 
and using the current and electric field as boundary conditions 
at the interface. The results of the simulation show that the 
holes pile up at the injecting layer side of the interface creat- 
ing a large electric field, which remains constant across the ac- 
cepting layer. In the range of applied voltage considered here, 
there is a low and constant concentration of holes throughput 
the accepting layer. The current in the accepting layer is sim- 
ply given by J = eiip int2 F int . 

We have modified the injection model to account for the 
hole occupation by summing the contribution to the current 
of all the occupied hole states, pi nt {(E), at the injecting layer 
side of the interface: 
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where N t is the total DOS, e is the electron charge, Vo is the 
attempt-to-jump frequency of the order of the phonon fre- 
quency, £j nt i is the energy of carriers at the injecting layer 
side of the interface, a is the nearest-neighbor distance, y is 
the inverse localization radius, and x is the hopping distance 
into the accepting layer. E m a is the energy of the first site on 
the accepting side, </>o is the barrier height, f[ nt is the interface 
electric field, g(E) =N t /V2na exp(-£' 2 /2a 2 ) is the Gaus- 
sian DOS in the accepting layer, and Bol(E) is the energy 
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The expression inside the square bracket in Eq. (1), without 
the term £i n ti, was derived by Arkhipov et al [11]. It calcu- 
lates the thermally assisted hopping rate for a hole at the inter- 
face of the injecting layer to a site over a distance, x, into the 
accepting layer. This is multiplied by the escape probability, 
Wesc, of the hole continuing into the bulk of the accepting layer 
(as opposed to hopping back down the barrier). The expres- 
sion from Ref. [11] assumes that there is a hole, at the HOMO 
level of the interface of the injecting layer, for every receiving 
site in the accepting layer. Our modification accounts for the 
occupancy and the energy distribution of holes at the injecting 
layer side of the interface. In order to do this, the expression 
in the square bracket is multiplied by the probability that an 
injecting site is occupied, Pma{E) /N t , and everything is inte- 
grated over Zijnti . The argument of Bo\(E) was modified from 
( <j)o + E int2 - eFmtx) in Ref . [1 1] to (0o - E mt i + E mt2 - eF^x) 
to account for the energy distribution in the injecting layer. 

Using the ATLAS simulation results for pi nt i, the quasi- 
Fermi level, Ep, and the hole occupation probability p mt \ can 
be obtained using 



Pintl=/ AEp m \(E) = 



N, 



2%a 



AE 



1 + exp ( 



E—Ep 
kT 



(3) 



where a is the width of the Gaussian DOS in the injecting 
layer. p- mt \ and the integrand of Eq. (1), J(E ml ), are plotted as 
a function of E^ nt in Fig. 1 for different values of Ep- Remark- 
ably, the current contributed by those carriers at Ep is almost 
negligible. The main contribution is due to the carriers at the 
energy of the maximum of the DOS, even though the Fermi 
level and the majority of the carriers lie below this level. This 
demonstrates that the assumption made in Ref. [12] that the 
carriers are injected from the Fermi level is incorrect. Mathe- 
matically, this peak is due to the exponentially falling factor of 
the Fermi function in p; nt i (E) in Eq. (3) canceling the expo- 
nentially rising factor in Bol(E) in Eq. (2), leaving the Gaus- 
sian in pi nt i(E) as the predominant energy dependent factor. 

The model described above assumes a constant energy level 
alignment at the organic heterojunction which is determined 
by imposing vacuum level alignment. However, the existence 
of an interface dipole would shift the vacuum level and change 
the alignment [18]. As discussed earlier in this letter, there 
is recent experimental and theoretical support for the exis- 
tence of an interface dipole in, at least, some organic systems. 
Considering this, we propose a model where the energy level 
alignment is determined by the changes in Fermi levels in the 
two organic layers. The result of this is that a change in car- 
rier concentrations will cause a change in the barrier height. 
The simplest assumption is to assume thermodynamic equi- 
librium across the interface which implies that the Fermi level 
is continuous at the interface. 
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FIG. 3: The field and voltage dependence of the variable barrier (j> v 
at 323 K (□), 288 K (O), 266 K (A), 219 K (y), and 172 K (O). 
Assuming thermodynamic equilibrium, 0,, is calculated using Eq. (3) 
with p mt \ and p m2 obtained by the steady-state simulation. The solid 
line is a fitting curve to an exponential decay function. 



This model is illustrated in Fig.2, which shows the energy 
levels in the case with an applied voltage and a current in the 
device. The energy level bending in the injecting layer is due 
to the holes piling up at the interface. Fig. 2(a) shows the 
case where the vacuum levels of the two layers are matched 
at the interface. In this case the Fermi level will be discontin- 
uous. The barrier height, which is the difference between the 
HOMO levels at the interface, is constant. Fig2(b) represents 
our model, where an interface dipole has formed which aligns 
the two Fermi levels. This will shift the vacuum levels and 
now the difference between the HOMO levels has shifted re- 
sulting in a smaller barrier height. Changing the current will 
change the hole concentrations, hence changing the Fermi lev- 
els. Evidently, this model features a variable barrier height 0,.. 

Using Eq. (3) and the simulated interface hole densities 
Pi ntl and p int2 of the ITO / MTDATA (320 nm) / NPB (390 
nm) / Ag heterojunction device, the variable barrier, r can be 
obtained by adjusting the energy alignment so that the Fermi 
levels E F match at the MTDATA/NPB interface. Our choice 
of the values of (Jmtdata and Onpb used in the calculation will 
be discussed later. Figure 3 shows the calculated 0,, as a func- 
tion of applied voltage at different temperatures. 0,, varies lit- 
tle with temperature, but it decreases exponentially when the 
applied voltage increases. This can be explained by the occu- 
pancy of holes in the interface DOSs. According to the sim- 
ulation results, with increasing voltage, p mt \ increases from 
10 17 cm~ 3 to 10 19 cm~ 3 and pi nt2 increases from 10 9 cm~ 3 to 
10 13 cm~ 3 . The amount of change of the Fermi levels depends 
on the changes in and p mt i and pi nt2 on the DOS (at the posi- 
tion of the Fermi levels). In this device, the Fermi level in the 
accepting layer changes more rapidly than that in the inject- 
ing layer, leading to a decreasing barrier as voltage increases. 
It can be deduced that, besides the energy difference of the 
HOMO levels, the energy profile is significant in controlling 
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FIG. 4: Measured J - F int characteristics of the ITO / MTDATA(320 
nm) / NPB(390 nm) / Ag heteroj unction device at 323 K (□), 288 
K (O), 266 K (A), 219 K (y), and 172 K (O). The correspond- 
ing solid symbols are the calculation results based on the modified 
injection model in Eq. (1) with the variable barriers, <j) v , in Fig. 3. 
The inset is the J vs \/T plot at Fj„ t = 0.5 M V/cm. Fitting the 
calculation results (solid line) with Eq. (1) to the experimental data 
(opened symbols) obtains (j) v = 0.45 ±0.1 eV, Omtdata = 95 meV 
and 0"npb = HOmeV. 

carrier injection across organic heterojunction. 

The inset in Fig. 4 illustrates the first step in our calcula- 
tion for the ITO / MTDATA / NPB / Ag heterojunction de- 
vice. An inverse localization radius y = 5 x 10 7 cm _1 has 
been used [19]. The nearest-neighbor distance a = 1 nm is 
taken as the size of a NPB molecule. By assuming only one 
net charge is carried by a molecule at a given time, the to- 
tal number of DOS N t = 1 x 10 21 cm~ 3 . For a fixed inter- 
face electric field (Fi nt = 0.5 MV/cm), the disorder parameters 
Omtdata an d <7npb> an d the attempt-to-jump frequency Vo, are 
varied to obtain the best fit to the temperature dependent data. 
Using the values v = 7.88 x 10 13 s~\ <j> v = 0.45 ± 0.1 eV, 
Omtdata = 95meV and Onpb = llOmeV, the theory agrees 
with the observed thermally activated current. The values for 
the disorder are slightly larger than those obtained from the 
TOF measurement. A broadening would be expected in the 
presence of an interface dipole [20]. The above values are 
used in Eq. (1), with the introduction of the variable barrier, 
<j> v , to calculate the J vs fi nt characteristics in Fig. 4. The 
calculated electric field and temperature dependent injection 
current densities are in excellent agreement with the experi- 
mental results. The good agreement between two sweeping 
directions of the applied voltage at all temperatures also sug- 
gests that the trap charging effect is negligible. It is worth- 
while noting that ignoring the contribution of the variable bar- 
rier, 0,„ will vastly underestimate the field dependence of the 
injection current across the heterojunction. 

In conclusion, an injection model for an organic heterojunc- 
tion which includes the distribution of carriers at the injecting 
side of the interface has been proposed. This model shows that 



the main contribution to the injection current is from carriers 
at the peak of the Gaussian density-of-state (DOS) rather than 
at the Fermi level. In addition, we have proposed a variable in- 
jection barrier height, 0,, which arises from the energy level re- 
alignment triggered by the relative Fermi level shifting caused 
by charge carrier accumulation at the interface. The model 
is verified by the excellent agreement with the experimental 
results of the current density voltage (J — V) characteristics 
of an organic heterojunction device for various temperatures 
over a wide range of electrical field. The results show that 
a dynamic change in carrier concentration and DOS energy 
alignment at the heterojunction plays a significant role in the 
performance of organic heterojunction devices. 
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